The Stochastic Dynamics of Rectangular and V-shaped Atomic Force Microscope 
Cantilevers in a Viscous Fluid and Near a Solid Boundary 
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Using a thermodynamic approach based upon the fluctuation-dissipation theorem we quantify the 
stochastic dynamics of rectangular and V-shaped microscale cantilevers immersed in a viscous fluid. 
We show that the stochastic cantilever dynamics as measured by the displacement of the cantilever 
tip or by the angle of the cantilever tip are different. We trace this difference to contributions 
from the higher modes of the cantilever. We flnd that contributions from the higher modes are 
signiflcant in the dynamics of the cantilever tip-angle. For the V-shaped cantilever the resulting 
flow field is three-dimensional and complex in contrast to what is found for a long and slender 
rectangular cantilever. Despite this complexity the stochastic dynamics can be predicted using a 
two-dimensional model with an appropriately chosen length scale. We also quantify the increased 
fluid dissipation that results as a V-shaped cantilever is brought near a solid planar boundary. 
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I. INTRODUCTION 

The stochastic dynamics of micron and nanoscale can- 
tilevers immersed in a viscous fluid are of broad scientific 
and technological interest [H, i]. Of particular impor- 
tance is the oscillating cantilever that is central to atomic 
force microscopy [3, Significant theoretical progress 
has been made using simplified models in the limit of long 
and thin rectangular cantilevers [1, H, 01- In this case, a 
two-dimensional approximation is appropriate (therefore 
neglecting effects due to the tip of the cantilever) and 
has yielded important insights. However, it is not cer- 
tain how well these approximations work for many sit- 
uations of direct experimental interest. For example, a 
commonly used cantilever in atomic force microscopy is 
V-shaped and a theoretical description of the dynamics 
of these cantilevers in fluid is not available. 

Furthermore, micron and nanoscale cantilevers are of- 
ten used in close proximity to a solid boundary either by 
necessity or out of experimental interest, ft is well known 
experimentally and theoretically that the presence of a 
solid boundary increases the fluid dissipation resulting 
in reduced quality factors and reduced resonant frequen- 
cies 0, i, H [13, im, [H, El, El- Again, theoretical de- 
scriptions are available in the limit of long and thin rect- 
angular cantilevers and it is uncertain if these approaches 
can be applied to these more complex geometries. 

In this paper we use a powerful thermodynamic ap- 
proach to quantify the stochastic dynamics of cantilevers 
due to Brownian motion for experimentally relevant ge- 
ometries for the precise conditions of experiment includ- 
ing the presence of a planar boundary. Our results are 
valid for the precise three-dimensional geometry of inter- 
est and include a complete description of the fluid-solid 
interactions. Using these results we are able to compare 
with available theory to yield further physical insights 



and to suggest simplified analytical approaches to de- 
scribe the cantilever dynamics for these complex situa- 
tions. 



II. THERMODYNAMIC APPROACH - 
FLUCTUATIONS FROM DISSIPATION 

The stochastic dynamics of micron and nanoscale can- 
tilevers driven by thermal or Brownian motion can be 
quantified using strictly deterministic calculations. This 
is accomplished using the fluctuation-dissipation theorem 
since the cantilever remains near thermodynamic equi- 
librium ^15; , d6j . We briefly review this approach for the 
case of determining the stochastic displacement of the 
cantilever tip and then extend it to the experimentally 
important case of determining the stochastic dynamics 
of the angle of the cantilever tip. 

The autocorrelation of equilibrium fluctuations in can- 
tilever displacement can be determined from the deter- 
ministic response of the cantilever to the removal of a 
step force from the tip of the cantilever (i.e. a transverse 
point force removed from the distal end of the cantilever) . 
If this force f{t) is given by 



fit) = 



Fq for i < 
for t > 0, 



(1) 



where t is time and Fq is the magnitude of the force, then 
the autocorrelation of the equilibrium fluctuations in the 
displacement of the cantilever tip is given directly by 



(ui(O)ui(i)) = keT 



Fn ' 



(2) 
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where is Boltzmann's constant, T is the tempera- 
ture, and is an equilibrium ensemble average. In our 
notation lower case letters represent stochastic variables 
{ui{t) is the stochastic displacement of the cantilever tip) 
and upper case letters represent deterministic variables 
{Ui{t) represents the deterministic ring down of the can- 
tilever tip due to the step force removal). The spectral 
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FIG. 1: Schematics of the two micron scale cantilever geome- 
tries considered (not drawn to scale). Panel (a), A rectan- 
gular cantilever with aspect ratios L/h = 98.5, w/h — 14.5, 
and L/w — 6.8. The cantilever is composed of silicon with 
density pc = 2329 kg/m^ and Youngs Modulus £ = 174 
GPa. Panel (b), A V-shaped cantilever with aspect ratios 
L/h = 233, w/h = 30, and L/w = 7.8. The total width 
between the two arms normalized by the width of a single 
arm is b/w — 10.36. The cantilever planform is an equilateral 
triangle with 9 — n/3. The cantilever is composed of silicon 
nitride with = 3100kg/m^ and E = 172GPa. The specific 
dimensions for the rectangular and V-shaped cantilever are 
given in Table |T] 



properties of the stochastic dynamics are given by the 
Fourier transform of the autocorrelation. 

The thermodynamic approach is valid for any conju- 
gate pair of variables 15]. For example, it is common 
in experiment to use optical techniques to measure the 
angle of the cantilever tip as a function of time [3]. It 
has also been proposed to use piezoresistive techniques 
to measure voltage as a function of time [l7[- The ther- 
modynamic approach remains valid for these situations 
by choosing the correct conjugate pair of variables. 

In this paper we also explore the stochastic dynamics 
of the angle of the cantilever tip. In this case, the angle 
of the cantilever tip is conjugate to a step point-torque 
applied to the cantilever tip. If this torque is given by 



r(t) 



To for i < 
for i > 0, 



(3) 



(1) 197 29 

(2) 140 15.6 



2 

0.6 



1.3 
0.1 



1.6 X 10^* 71 
8.9 X 10"^° 38 



TABLE I: Summary of the cantilever geometries and material 
properties. (1) The rectangular cantilever. (2) The V-shaped 
cantilever used is the commercially available Veeco MLCT 
Type E microlever that is used in AFM [l^ . The geometry 
is given by the cantilever length L, width w, and height h. 
For the V-shaped cantilever the total length between the two 
arms at the base is 6 = 161.64^m. The cantilever spring con- 
stant k, torsional spring constant kt, and resonant frequency 
in vacuum fo are determined using finite element numerical 
simulations. The cantilevers are immersed in water with den- 
sity pi = 997 kg/m^ and dynamic viscosity ri — 8.59 x 10~* 
kg/m-s. 



tip- angle 9{t) is given by 



ei(t) 



To 



(4) 



Here Qi(t) represents the deterministic ring down, as 
measured by the tip-angle, resulting from the removal 
of a step point-torque. Again, the Fourier transform of 
the autocorrelation yields the noise spectrum. 

A powerful aspect of this approach is that it is possible 
to use deterministic numerical simulations to determine 
Ui{t) and Qi{t) for the precise cantilever geometries and 
conditions of experiment. This includes the full three- 
dimensionality of the dynamics which are not accounted 
for in available theoretical descriptions. The numerical 
results can be used to guide the development of more 
accurate theoretical models. 



III. THE STOCHASTIC DYNAMICS OF 
CANTILEVER TIP-DEFLECTION AND 
TIP-ANGLE 

The stochastic dynamics of the cantilever tip- 
displacement ui{t) and that of the tip-angle Oi{t) yield 
interesting differences. Using the thermodynamic ap- 
proach, insight into these differences can be gained by 
performing a mode expansion of the cantilever using the 
initial deflection required by the deterministic calcula- 
tion. The two cases of a tip-force and a tip-torque re- 
sult in a significant difference in the mode expansion co- 
efficients which can be directly related to the resulting 
stochastic dynamics. 

For small deflections the dynamics of a cantilever 
with a non-varying cross section are given by the Euler- 
BcrnouUi beam equation. 



at2 



(5) 



where tq is the magnitude of the step torque, then the 
autocorrelation of equilibrium fluctuations in cantilever 



where U{x,t) is the transverse beam deflection, fj, is the 
mass per unit length, E is Young's modulus, and / is the 



moment of inertia [l9l|. For the case of a cantilever where 
a step force has been appHed to the tip at some time in 
the distant past the steady deflection of the cantilever at 
i = is given by 



U{x) 



'2EI 



Lx' 

3 



(6) 



where L is the length of the cantilever and the appropri- 
ate boundary conditions are C/(0) = [/'(O) = U"{L) = 
and U"'{L) — —Fq/EI. The prime denotes differentia- 
tion with respect to x. 

Similarly, the deflection of the same cantilever beam 
due to the application of a point-torque at the cantilever- 
tip is quadratic in axial distance and is given by 



U{x) 



To 2 

2EI ' 



(7) 



where the appropriate boundary conditions are U{0) = 
U'{0) = U"'{L) = and C7"(L) = tq/EI. The angle 
of the cantilever measured relative to the horizontal or 
undisplaced cantilever is then given by tan 8 = U'{x). 
The mode shapes for a cantilevered beam are given by 

^nix) — — (cos + cosh kX) (cos Kx — cosh Kx) 

— (sin kL — sinh kL) (sin kx — sinh kx) , (8) 

where n is the mode number, and the characteristic fre- 
quencies are given by = us^^/EI. The mode numbers 
K are solutions to 1 + cos kL cosh kL — [31 • The initial 
cantilever displacement given by Eqs. ([6]) and ^ can be 
expanded into the beam modes 



Uix) 



a„$„(x). 



(9) 



with mode coeflicients a„. The total energy E\^ of the 
deflected beam is given by 



Ey, 



EI 



U"{xfdx, 



(10) 



which is entirely composed of bending energy. The frac- 
tion of the total bending energy contained in an individ- 
ual mode is given by 



EI 
2Eb 



{an'i>'^{x))- dx. 



(11) 



The coefficients 6„ for the rectangular cantilever of Ta- 
bleUare shown in Table [Til For the case of a force applied 
to the cantilever tip, 97% of the total bending energy 
is contained in the fundamental mode and the energy 
contained in the higher modes decays rapidly with less 
than 1% of the energy contained in mode three. When 
a point-torque is applied to the same beam it is clear 
that a significant portion of the bending energy is spread 
over the higher modes. Only 61% of the energy is con- 
tained in the fundamental mode and the decay in energy 



n 


On (tip-iorce) 


On (tip-torque) 


1 


0.97068 


0.61308 


2 


0.02472 


0.18830 


3 


0.00315 


0.06473 


4 


0.00082 


0.03309 


5 


0.00030 


0.02669 



TABLE II: The fraction of the total energy E^y contained in 
the first five beam modes given by the coeflicients 6„. The 
tip-force results are for a rectangular beam that has been 
deflected by the application of a point force to the cantilever 
tip. The tip-torque results are for a rectangular beam that 
has been deflected by the application of a point torque to 
the cantilever tip. The coefficients clearly show that the tip- 
torque case has significantly more energy contained in the 
higher modes. 



1/2, 



(1) 

(2) 



5.6 

20 



5.0 X 10 
7.0 X 10 



TABLE III: The magnitude of stochastic fiuctuations in tip- 
defiection and in tip-angle for the rectangular (1) and V- 
shaped (2) cantilevers. These values were obtained from nu- 
merical simulations simulations of the beams in vacuum. 



with mode number is more gradual. The fifth mode for 
the tip-torque case contains more energy than the second 
mode for the tip-force case. Although we have only dis- 
cussed a mode expansion for the rectangular cantilever, 
the V-shaped cantilever will exhibit similar trends since 
the transverse mode shapes are similar to that of a rect- 
angular beam [20t . 

The variation in the energy distribution among the 
modes required to describe the initial deflection of the 
cantilever can be immediately connected to the resulting 
stochastic dynamics. For the deterministic calculations 
the initial displacement can be arbitrarily set to a small 
value. In this limit the modes of the cantilever beam 
are not coupled through the fluid dynamics. As a result, 
the stochastic dynamics of each mode can be treated as 
the ring down of that mode from the initial deflection. 
This indicates that the more energy that is distributed 
amongst the higher modes initially the more significant 
the ring down and, using the fluctuation-dissipation the- 
orem, the more significant the stochastic dynamics. 

The mode expansion clearly shows that the tip-torque 
case has more energy in the higher modes. This suggests 
that stochastic measurements of the cantilever tip-angle 
will have a stronger signature from the higher modes than 
measurements of cantilever tip-displacements. Using fi- 
nite element simulations for the precise geometries of in- 
terest we quantitatively explore these predictions. 
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IV. THE STOCHASTIC DYNAMICS OF A 
RECTANGULAR CANTILEVER 

We have performed deterministic numerical simula- 
tions of the three-dimensional, time dependent, fluid- 
solid interaction problem to quantify the stochastic dy- 
namics of a rectangular cantilever immersed in water 
using the thermodynamic approach discussed in Sec- 
tionini The deterministic numerical simulations are done 
using a flnite element approach that is described else- 
where [211^22,]. 

The stochastic fluctuations in cantilever tip- 
displacement for a rectangular cantilever in water 
have been described elsewhere [l^ [l^ [H, In the 
following we compare these results with the stochastic 
dynamics as determined by the fluctuations of the 
cantilever tip-angle. The geometry of the the specific 
micron scale cantilever we explore is given in Table HI 

As discussed in Section |lT] the autocorrelations in equi- 
librium fluctuations follow immediately from the ring 
down of the cantilever due to the removal of a step 
force (to yield (ui(O)wi(t))) or step point-torque (to yield 
(6*1 (0)01 (t))). The autocorrelations of the rectangular 
cantilever are shown in Fig. [51 The magnitude of the 
noise is quantifled by the root mean squared tip-angle 
and deflection which is listed in Table [TTTl 

A comparison of the autocorrelations yields some in- 
teresting features. At short times {9i{0)9i{t)) shows the 
presence of higher harmonic contributions. This is shown 
more clearly in the inset of Fig. [51 This further sug- 
gests that the angle autocorrelations are more sensitive 
to higher mode dynamics as discussed in Section [TTTl 




^0 0.05 0.1 0.15 0.2 

t(ms) 



FIG. 2: The normalized autocorrelation of the rectangular 
cantilever for tip-deflection (solid) and tip-angle (dashed). 
(Inset) A detailed view of the autocorrelation at short time 
differences to illustrate the influence of higher modes in the 
tip-angle measurements. 

The Fourier transform of the autocorrelations yield the 



noise spectra shown in Fig. [3l In our notation the sub- 
script of G indicates the variable over which the noise 
spectrum is measured: Gg is the noise spectrum for tip- 
angle and Gu is the noise spectrum for tip-displacement. 
The equipartition theorem of energy yields, 

- Geiu;)dcj = -j- (13) 

where k and kt are the transverse and torsional spring 
constants, respectively. The curves in Fig. [31 are normal- 
ized using the equipartition result to have a total area 
of unity. Using this normalization the area under a peak 
is an indication of the amount of energy contained in 
a particular mode. Figure [31 shows only the flrst two 
modes, although the numerical simulations include all of 
the modes (within the numerical resolution of the finite 
element simulation). The energy distribution across the 
first two modes shows the significance of the second mode 
for the tip-angle dynamics. 
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FIG. 3: The noise spectra of stochastic fluctuations in can- 
tilever tip-angle (dashed) and tip-deflection (solid) for the 
rectangular cantilever. The curves are normalized to have 
the same area, however only the first two modes are shown. 

Using a simple harmonic oscillator approximation it is 
straight forward to compute the peak frequency ujf and 
quality Q for the cantilever in fluid. Using a single mode 
approximation yields the values shown in Table IIVI As 
expected there is a significant reduction in the cantilever 
frequency when compared with the resonant frequency 
in vacuum ujq and the quality factor is quite low because 
of the strong fluid dissipation. The values of tj/ and Q 
for tip-angle and tip-dispacement are nearly equal. This 
is expected since the displacements and angles are very 
small, resulting in negligible coupling between the modes. 
Any differences in and Q can be attributed to using 
a single mode approximation. 
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It is useful to compare these results with the commonly 
used approximation of an oscillating, infinitely long cylin- 
der with radius w/2 [1, H, The cantilever used here 
has an aspect ratio of L/w ~ 7 and the infinite cylinder 
theory is quite good at predicting of and Q. 

LOf/iOO Q 



(1) 0.35 3.34 

(2) 0.36 3.26 

TABLE IV: The peak frequency and quality factor of the 
fundamental mode of the rectangular cantilever determined 
by finite element simulations using the thermodynamic ap- 
proach. (1) is computed using the cantilever tip-displacement 
due to the removal of a step force. (2) is computed using the 
cantilever tip-angle due to the removal of a point-torque. The 
frequency result is normalized by the resonant frequency in 
vacuum uq. Using the infinite cylinder approximation with 
a radius of to/2 the analytical predictions are Q — 3.24 and 
LUf/ujQ = 0.34. 



V. THE STOCHASTIC DYNAMICS OF A 
V-SHAPED CANTILEVER 

We now explore the stochastic dynamics of a V-shaped 
cantilever in fluid. An integral component of any theo- 
retical model is an analytical description of the result- 
ing fluid flow flcld caused by the oscillating cantilever. 
The deterministic finite element simulations that we per- 
formed yield a quantitative picture of the resulting fiuid 
dynamics. Exploring the flow fields further yields insight 
into the dominant features that contribute to the can- 
tilever dynamics. 

As discussed earlier, for long and slender rectangular 
cantilevers the flow field is often approximated by that 
of a cylinder of diameter w undergoing transverse oscil- 
lations. This approach assumes that the fiuid fiow is es- 
sentially two-dimensional in the y — z plane and neglects 
any fiow over the tip of the cantilever. Figure H] (top) il- 
lustrates this tip flow for the rectangular cantilever using 
vectors of the fluid velocity in the x — y plane at 2 = 0. 
The flgure is a close-up view near the tip of the can- 
tilever. It is evident that the flow over the rectangular 
cantilever is nearly uniform in the axial direction leading 
up to the tip. However, near the tip there is a significant 
tip flow that decays rapidly in the axial direction away 
from the tip. The increasing significance of the tip fiow 
as the cantilever geometry becomes shorter (for exam- 
ple, by simply decreasing L) is not certain and remains 
an interesting open question. However, for the geometry 
used here it is clear that this tip-fiow is negligible based 
upon the accuracy of the analytical predictions using the 
two-dimensional model. 

Figure S] (bottom) illustrates the tip flow for the V- 
shaped cantilever, again by showing velocity vectors in 
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FIG. 4: The fluid flow near the tip of the cantilever as il- 
lustrated by the velocity vector field calculated from finite 
element numerical simulations. A cross section of the x — y 
plane at z = is shown (see Fig. [T]) that is a close-up view 
of the tip-region. The shaded region indicates the cantilever 
(because of the small defiections used in the simulations that 
cantilever does not appear to be defiected). (top) The fiow 
field near the tip of the rectangular cantilever. This fiow field 
is at t=6/is and the magnitude of the largest velocity vector 
shown is -0.3 nm/s. (bottom) The fiow field near the tip of 
the V-shaped cantilever. This fiow field is at t=7.2/^s and the 
magnitude of the largest velocity vector shown is -26 nm/s. 
The shaded region indicates the tip region where the two sin- 
gle arms have merged. The open region to the left is where 
the two single arms have separated revealing the open region 
in the interior of the V-shaped cantilever. 



the x — y plane at z = 0. The shaded region indicates the 
part of the cantilever where the two arms have merged. 
To the right of the shaded region illustrates flow off the 
tip and to the left indicates flow that circulates back in 
between the two individual arms. 

In order to illustrate the three-dimensional nature of 
this flow, the flow field in the y — z plane is shown at two 
axial locations in Fig.[5l Figure[5l^top) is at axial location 
X = 77 fj-m. The two shaded regions indicate the two arms 
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of the cantilever. Each arm is generating a flow with a 
viscous boundary layer (Stokes layer) as expected from 
previous work on rectangular cantilevers. However, the 
Stokes layers interact in a complicated manner near the 
center. It is expected that as one goes from the base of 
the cantilever to the tip that these fluid structures would 
transition from non-interacting to strongly-interacting. 

Figure [S^bottom) illustrates the flow fleld at axial lo- 
cation X = 108.8/im, the axial location at which the two 
arms of the cantilever merge to form the tip region. The 
length of the shaded region is therefore 36/im or twice 
that of a single arm shown in Fig. [Iftop). For this tip 
region the flow field is similar to what would be expected 
of a single rectangular cantilever of this width. 

Overall, it is clear that the fluid flow field is more com- 
plex for the V-shaped cantilever than for the long and 
slender rectangular beam. For the V-shaped cantilever 
the flow is three-dimensional near the tip region where 
the two arms join together. 

Central to the flow field dynamics are the interactions 
of the two Stokes layers caused by the oscillating can- 
tilever arms. The thickness of these Stokes layers are 
expected to scale with the frequency of oscillation as 
5s/a ^ Ruj~^^^ where a is the half-width of the can- 
tilever and — wa^/v is a frequency based Reynolds 
number (often called the frequency parameter). For the 
relevant case of a cylinder of radius a oscillating at fre- 
quency Lo the solution to the unsteady Stokes equations 
yields a distance of approximately 5^^ to capture 99% of 
the fluid velocity in the viscous boundary layer [25| . For 
a single arm of the V-shaped cantilever this distance is 
nearly 10/im. In comparison, the total distance between 
the two arms at the base is 125/im. This separation is 
large enough such that the two Stokes layers have neg- 
ligible interactions near the base. However, as the arms 
approach one another with axial distance the Stokes lay- 
ers overlap and eventually merge at the tip. 

Despite the complicated interactions of the three- 
dimensional flow caused by the cantilever tip and the 
axial merging of the two Stokes layers, the V-shaped can- 
tilever behaves as a damped simple harmonic oscillator. 
The autocorrelations in tip-angle and tip-displacement 
that are found using full finite element numerical sim- 
ulations are shown in Fig. [Sj It is again clear that the 
tip-angle dynamics have significant contributions from 
the higher modes, see the inset of Fig. [51 The area nor- 
malized noise spectra are shown in Fig. [71 

Using a simple harmonic oscillator analogy a peak fre- 
quency and a quality factor can be determined from the 
first mode in the noise spectra of Fig. [71 These values 
are given in the first two rows of Table IVII The qual- 
ity of the cantilever is Q « 2 and the peak frequency is 
reduced significantly, ujf/tijQ k, 0.2. compared to the res- 
onant frequency in the absence of a surrounding viscous 
fiuid. 

It is insightful and of practical use to determine the 
geometry of the equivalent rectangular beam that would 
yield the precise values of k, ujf, and Q calculated for 
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FIG. 5: The fluid velocity vector field at two Eixial positions 
along the V-shaped cantilever calculated from deterministic 
finite element numerical simulations. Cross sections of the 
y—z plane are shown (see Fig.[T}, the entire simulation domain 
is not shown and the shaded region indicates the cantilever. 
Both images are taken at t=7.2/is and the maximum velocity 
vector shown is -26 nm/s. (top) The y — z plane ai x = 
77/im. The skewed width of a single arm of the cantilever 
in this cross-section is 18/im. The distance separating the 
two cantilever arms is 36^m. (bottom) The y — z plane at 
X = 108.8/im. This is the point at which the two single arms 
join to make a continuous cross-section of width 2w. 



the V-shaped cantilever from full finite-element numeri- 
cal simulations. For the rectangular beam the equations 
are well known (c.f. Ref. (l6j ) and yield a unique value 
of length L', width w' , and height h' as shown below, 



3^;/ Ew'h'^ 



Q 



L'3 

TrifOJf ^ TvpfW 



i£^-^r'(u;',c^/) 



7/ 



(14) 
(15) 



where the peak frequency is determined from the maxi- 
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mass, 7/ is the fluid damping, T is the hydrodynamic 
function for an infinite cyhnder, F' is the real part of F, 
and F" is the imaginary part of F. Equations (|14p - (fT5|) 
can be solved to yield values for the unknown geometry 
of the equivalent rectangular beam i', w' , and h' which 
are given in Table |Vl The equivalent beam is shorter, 
thinner, and wider than the V-shaped cantilever. Impor- 
tantly, the width of the equivalent beam is nearly twice 
that of a single arm of the V-shaped cantilever. 



L'/L 



O.g 



w' /w 



1.9 



h/h' 



0.8 
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FIG. 6: The normalized autocorrelation of equilibrium fluc- 
tuations in the tip-deflection {ui{0)ui{t)} (solid lined) and 
in tip-angle {6i{0)9i{t)) (dashed-line) for the V-shaped can- 
tilever. The inset shows a close-up of the dynamics for short 
time differences to illustrate the influence of the higher modes 
in the tip-angle measurements. 




FIG. 7: The noise spectra for the V-shaped cantilever as de- 
termined from the tip-displacement Gx (solid line) and from 
tip-angle Ge (dashed line). The curves are normalized to have 
an area of unity, with only the first two modes shown. 



TABLE V: The geometry of the equivalent rectangular beam 
that yields the exact values of k, uif, and Q for the V-shaped 
cantilever that have been determined from full finite-element 
numerical simulations. The length, width, and height of the 
equivalent beam {L' ,w' ,h') are calculated using Eqs. (|14|l - 
(|15p and are normalized by the values of {L,w,h) for the 
V-shaped cantilever given in Table |T] 

These results sug gest that the parallel beam approx- 
imation (PBA) [11113, m, m commonly used to de- 
termine the spring constant for a V-shaped cantilever 
may also provide a useful geometry for determining the 
dynamics of V-shaped cantilevers in fiuid. In this ap- 
proximation the V-shaped cantilever is replaced by an 
equivalent rectangular beam of length L, width 2w, and 
height h to yield a simple analytical expression for the 
spring constant. This has been shown to be quite suc- 
cessful for V-shaped cantilevers that have arms that are 
not significantly skewed. The results of using the geom- 
etry of this approximation to determine lo f and Q from 
the two-dimensional cylinder approximation are shown 
on the third row of Table IVTl It is clear that this is quite 
accurate. It is expected that these results will remain 
useful for cantilever geometries that do not deviate sig- 
nificantly from that of an equilateral triangle as studied 
here. An exploration of the breakdown of this approx- 
imation is possible using the methods described but is 
beyond the scope of the current efforts. 



VI. QUANTIFYING THE INCREASED 
DISSIPATION DUE TO A PLANAR BOUNDARY 



mum of the noise spectrum, 
4fcsT 1 



Gu — 



k ujQ 



(16) 



To[bT"{Ro[;;) 



[(1 - (^2(1 + ToT'iRoCuW + (^2ToF"(i?oc^))2] ' 

In the above equations uj — lo/luq is the normalized fre- 
quency, a — 0.234 is a constant factor to determine an 
equivalent lumped mass for a rectangular beam, is 
the equivalent mass of the cantilever plus the added fluid 



In practice, the cantilever is never placed in an un- 
bounded fluid and the influence of nearby boundaries 
must be accounted for to provide a complete description 
of the dynamics. In many cases the cantilever is purpose- 
fully brought near a surface out of experimental interest 
in order to probe some interaction with the cantilever or 
to probe the surface itself. To specify our discussion we 
will consider the situation depicted in Fig. [8] showing a 
cantilever a distance s from a planar boundary. In the 
following we study the case where the cantilever exhibits 
flexural oscillations in the direction perpendicular to the 
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1 




(1) 


0.21 


1.98 


(2) 


0.22 


2.04 


(L,2w,h) 


0.19 


1.98 



TABLE VI: The peak frequency and quality factor of the fun- 
damental mode of the V-shaped cantilever determined by fi- 
nite element simulations using the thermodynamic approach. 
(1) is computed using the cantilever tip-displacement due to 
the removal of a step force. (2) is computed using the can- 
tilever tip-angle due to the removal of a point-torque. The 
third line represents theoretical predictions using the geom- 
etry of an equivalent rectangular beam given by {L,2w,h). 
The frequency result is normalized by the resonant frequency 
in vacuum too- 

boundary. However, we would like to emphasize that 
our approach is general and can be used to explore ar- 
bitrary cantilever orientations and oscillation directions 
if desired. The fluid is assumed to be unbounded in all 
other directions. It is well known that the presence of the 
boun dary will influence the dynamics of the cantilever 
[1, [§, ITsf . The result is a reduction in the resonant fre- 
quency and quality factor. This has been described the- 
oretically for the case of a long and thin cantilever of 
simple geometry where the flui d dy namics have been as- 
sumed two-dimensional 0, HO; [HI • 




FIG. 8: A schematic of a cantilever a distance s away from 
a solid planar surface (not drawn to scale). The cantilever 
undergoes flexural oscillations perpendicular to the surface. 

In the following we use the thermodynamic approach 
with flnite element numerical simulations to quantify the 
dynamics of the V-shaped cantilever as a function of its 
separation from a boundary. We have performed 8 simu- 
lations over a range of separations from 10 to 60/im using 
both the tip-deflection and tip-angle formulations. The 
noise spectra for these simulations are shown in Fig. [9l 
Using the insights from our simulations of the V-shaped 
cantilever in an unbounded fluid we expect the relevant 
length scale for the fluid dynamics to be twice the width 
of a single arm, 2w. Using the peak frequency of the 
V-shaped cantilever in unbounded fluid yields a Stokes 
length Ss — 4.14/im. Scaling the separation by the Stokes 
length yields. 2.5 < s/Sg < 15 which covers the range 



from what is expected to be a strong influence of the 
wall to a negligible influence. Figure [9] clearly shows a 
reduction in the peak frequency and a broadening of the 
peak as the cantilever is brought closer to the boundary. 
In fact, for the smaller separations the peak is quite broad 
and the trend suggests that eventually the peak will be- 
come annihilated as the cantilever is brought closer to 
the boundary. 




0.1 0.2 , 0.3 0.4 0.5 

CO/ COq 



FIG. 9: Panel (a) The noise spectra Gu of stochastic fluc- 
tuations in cantilever tip-deflection for separations s — 
10, 12, 15, 20, 40/im. Panel (b) the noise spectra Ge of stochas- 
tic fluctuations in cantilever tip-angle for separations s = 
15, 25, 60/im. The spectra have been normalized by the max- 
imum value of Gu or Gq. The smallest and largest values of 
separation are labeled with all other values appearing sequen- 
tially. 

Using the noise spectra we compute a peak frequency 
and a quality factor for the fundamental mode as a func- 
tion of separation from the boundary, which are plot- 
ted in Fig. [ini The horizontal dashed line represents 
the value of the peak frequency and quality factor in the 
absence of bounding surfaces using the two-dimensional 
inflnite cylinder approximation [l6| where the cylinder 
width has been chosen to be 2w. It is clear from the 
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results that for separations greater than s/5s ^ 7 the 
V-shaped cantilever is not significantly affected by the 
presence of the boundary. However, as the separation 
decreases below this value the peak frequency and qual- 
ity factor decrease rapidly. 

The triangles in Fig. [10] represent the theoretical 
predictions of Green and Sader [111, |T3| using a two- 
dimensional approximation for a beam of uniform cross- 
section that accounts for the presence of the boundary. 
We have used a width of 2w in computing these theo- 
retical predictions for comparison with our numerical re- 
sults. Despite the complex and three-dimensional nature 
of the flow field the theory is able to accurately predict 
the quality factor over the range of separations explored. 
The frequency of the peak for the V-shaped cantilever 
shows some deviation from these predictions. 

In general, an increase in the period of oscillation for 
a submerged object can be attributed to the mass of 
fluid entrained by the object [13] ■ The lower peak fre- 
quency calculated for the V-shaped cantilever using a 
two-dimensional solution indicates an over-prediction of 
the mass loading. This can be attributed to the three- 
dimensional flow around the tip being neglected for this 
approach. It is reasonable to expect the cantilever tip to 
carry a smaller amount of fluid than a section of the beam 
body moving with the same velocity, see Fig.|4l The qual- 
ity factor relates to the ratio of the mass loading and the 
viscous dissipation and is less sensitive to deviations in- 
curred from the two-dimensional approximation. Despite 
neglecting three-dimensional flow around the cantilever 
tip, the two-dimensional model for the fluid flow around 
the V-shaped cantilever gives an accurate prediction of 
the peak frequency and quality factor. 





VII. CONCLUSIONS 

We have shown that the thermodynamic approach is a 
versatile and powerful method for predicting the stochas- 
tic dynamics of cantilevers in fluid for the precise con- 
ditions of experiment including complex geometries and 
the presence of nearby boundaries. Available analyti- 
cal predictions are for idealized situations including sim- 
ple geometries where the three-dimensional flow near the 
cantilever tip has been neglected. Although this has pro- 
vided signiflcant insight, many situations of experimen- 
tal interest are more complicated. It is often required to 
have a quantitative base-line understanding of the can- 
tilever dynamics for the precise conditions of experiment 
in order to make and interpret measurements in novel 
situations and in the presence of other phenomena of in- 
terest. 

We emphasize that by using the fluctuation-dissipation 
theorem a single deterministic calculation is sufficient to 
predict the stochastic behavior for all frequencies. Fur- 
thermore, the deterministic calculation is computation- 
ally inexpensive and does not require special computing 
resources. 



FIG. 10: The variation of the peak frequency (panel (a)) and 
quality (panel (b)) of the fundamental mode of the V-shaped 
cantilever in fluid as a function of separation from a nearby 
wall. Results calculated using tip-deflection are circles, results 
using tip-angle are squares, and theoretical predictions using 
the results of Ref. [Ij] are triangles. The peak frequency 
and quality factor of the fundamental mode in an unbounded 
fluid are u)/uj[ « 0.19 and Q ~ 2 and are represented by the 
horizontal dashed line. The distance s is normalized by the 
Stokes length 5s where a = w to yield 5s = 4.14/im. 

The thermodynamic approach is general in that it can 
be used to compute the stochastic dynamics of any conju- 
gate pair of variables. We have shown that the stochastic 
dynamics that are measured depend upon the choice of 
measurement. This could be exploited in future exper- 
iments, for example, to minimize or maximize the sig- 
nificance of the higher mode dynamics by choosing to 
measure tip-deflection or tip-angle, respectively. 

Our results also suggest that despite the complicated 
three-dimensional nature of the flow field around a V- 
shaped cantilever, analytical predictions based upon a 
two-dimensional description are surprisingly accurate if 
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the appropriate length scales are used. We anticipate 
that these findings will be of immediate use as the atomic 
force microscope continues to find further use in liquid 
environments. 
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